Angle Optimization of Graphs Embedded in the Plane 



Sergey Bereg* Timothy Rozario^ 

Abstract: In this paper we study problems of drawing graphs in the plane using edge length constraints 
and angle optimization. Specifically we consider the problem of maximizing the minimum angle, the MMA 
problem. We solve the MMA problem using a spring-embedding approach where two forces are applied to the 
vertices of the graph: a force optimizing edge lengths and a force optimizing angles. We solve analytically 
the problem of computing an optimal displacement of a graph vertex optimizing the angles between edges 
incident to it if the degree of the vertex is at most three. We also apply a numerical approach for computing 
the forces applied to vertices of higher degree. We implemented our algorithm in Java and present drawings 
of some graphs. 

1 Introduction 

Angular resolution is one of the aesthetic criteria measuring the quality of graph drawings in terms of human 
comprehension. The angular resolution of a straight-line drawing in the plane is the minimum angle between 
any two incident edges. The study of graph drawing with angular resolution started by Formann et aZ. jlOj 
in 1990. They introduced the angular resolution of a graph as the supremum angular resolution over all 
straight-line drawings of the graph. The problem of computing the angular resolution of a graph is NP-hard 
(even the problem of deciding if R = n/2 for graphs with vertex degrees at most four is NP-hard). 

The main focus in the early investigation [101 115] was on bounding the angular resolution of a graph 
in terms of the maximum vertex degree d. The obvious upper bound is R(d) < 2n/d. A lower bound 
R(d) > £1(1/ d) has been proved [10] for many graphs including planar graphs and complete graphs. A lower 
bound R(d) > £l(l/d 2 ) holds for all graphs [10]. If we insist on planar straight-line drawings, then R(d) > a d 
for some constant < a < 1 [15] , 

There was also a study of the optimization problem where the angular resolution of a given graph is 
maximized. Carlson and Eppstein [2] considered tree drawings such that all faces form (infinite) convex 
polygons. They showed that the optimal angular resolution can be computed in linear time and even 
the lengths of the tree edges may be chosen arbitrarily. In the recent paper, Eppstein and Wortman [9] 
considered graph drawing in the plane with faces drawn as centrally symmetric convex polygons. They 
found a polynomial time algorithm for computing a drawing maximizing the angular resolution. 

Another direction in graph drawing with the angle resolution is Lombardi drawing where the graph edges 
are represented as circular arcs instead of straight-line segments [3 El [7] . Relaxing the condition of straight- 
line segments allows to achieve perfect angular resolution where the edges are equiangularly spaced around 
each vertex. The classes of graphs admitting Lombardi drawings are presented in [7] (and the algorithms 
for finding these drawings). Duncan et al. [B] found that unrooted trees can drawn with perfect angular 
resolution and polynomial area. 

In this paper we also study the problem of maximizing the angular resolution using the force-directed 
approach. This idea is not new and some algorithms for optimizing angular resolution using the force-directed 
or spring-embedding approach [2] [5] . The common feature of these approaches is that the forces are directed 
from one vertex to another vertex. Our approach is different in the sense that we want to optimize the 
dislocation of a vertex. It can be viewed as a perturbation problem where a graph is embedded in the plane 
with straight-line edges and we want to find a new position of a given vertex by moving it at distance at 
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most r > (for a given r) and optimizing the incident angles. The perturbation problem is interesting in its 
own right. It is also interesting that the problem is related to the famous Format problem (which appears 
as a special case when the vertex degree is three). 

Another reason why we introduce and study the perturbation problem is that the parameter r allows to 
control the strength of the angular resolution force applied to vertices. This works well if we use more than 
one force. For example, we applied it to the spring embedding with length constraints. 

The paper is organized as follows. In Section [2] we briefly describe force-directed graph drawing and 
introduce Max-Min- Angle problem. In Section [3] we recall the classical Fermat problem that appears as the 
special case of the MAA problem. In Section [4] we provide a solution of the MM A problem for vertex degree 
two. Optimal solution for degree three vertices is provided in Section [5] The algorithm and its performance 
is discussed in Section [6] Finally we conclude in Section [7] 

2 Spring Embedder and Problem Statement 

Force-directed graph drawing is a popular technique and there is a growing literature on force-directed 
drawing algorithms, see the recent survey by Kobourov |14j . Eades [8] introduced a mechanical model for 
graph drawing. To achieve aesthetically pleasing layouts and capture the edge length constraints he applied 
attracting/repelling force between two vertices if the distance between them is less/greater than the desired 
length. He found that the Hookes Law (linear) springs are too strong when the vertices are far apart but 
the logarithmic force solves this problem. As initial embedding of the graph the algorithm places its vertices 
at random locations. The algorithm stops after a sufficient number of iterations. If a state of equilibrium is 
reached i.e. all forces are zero, then the graph embedding reaches the desired positioning in the plane and 
remains static. An example of such drawing is shown in Figure [T] 



Figure 1: A drawing of Algae graph 16J produced by the spring embedder pQ. 

Fruchterman and Reingold lij added a condition of "even vertex distribution" which is modeled by 
attractive forces between adjacent vertices and repulsive forces between all pairs of vertices. This increases 
the number of forces {\V\ 2 repulsive forces for a graph G = (V,E)) and slows down the algorithm. 

The algorithms by Eades [5] and Fruchterman and Reingold [TT] are just two examples of force-directed 
graph drawing. There are many spring embedders nowadays [14] . We consider a new force-directed approach 
for angular resolution where the vertex displacement is optimized locally. 

2.1 Problem Statement 

We consider one step of the spring embedder where the vertices of a graph G are embedded in the plane and 
they are allowed to move slightly. The spring embedder may take into account several forces, for example the 
forces aiming to achieve the desired lengths of the edges of G. We introduce a force that aims to optimize, 
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for all vertices v, the angles between edges incident to v and embedded in the plane. Ideally, we want the 
edges to spread evenly around v. This may not be possible due to other constraints of the drawing (edge 
lengths, for example). Our goal is to maximize, for every v, the smallest angle between two edges incident 
to v. For this, we formulate our problem. 

Max-Min- Angle (MMA) Problem. Let G be a graph embedded in the plane. Let P be the 
position of a vertex v of G and let A\, A2, . . . , Ak be the positions of the adjacent vertices of P 
in G. Let r > be a radius. Find a point P* (the best position to move P within distance r) 
such that I PP* I < r and the smallest angle ZAiP* Aj is maximized. 

Let Z be the circle centered at P and radius r. We will solve the MMA problem in Sections [4] and [5] by 
considering degree k of vertex v. 

3 Fermat problem 

When the degree of v is three, the MMA problem is related to the classical Fermat problem: Given triangle 
ABC in the plane, find a point F such that the total distance from the three vertices of the triangle to F is 
the minimum possible. The solution of the Fermat problem (called Fermat point or Torricelli point) depends 
on triangle ABC. 



B' 




(a) ( b ) 

Figure 2: The Fermat point F in (a) Case 1 and (b) Case 2. (a) Three angles ZAFB, ZBFC, and ZCFA 
are equal, (b) The Fermat point is at C. 

Case 1: All angles of triangle ABC are smaller than 120°. 
Construct three new regular triangles ABC' , AB'C, and A'BC out of the three sides of triangle ABC. Then 
the point F is the intersection of three lines AA',BB\ and CC, see Figure [2^a) . In this case the Fermat 
point F is coincident with the first isogonic center of the triangle (18] , where the angles subtended by the 
sides of the triangle are all equal, i.e. ZAFB, ZBFC, ZCFA = 120°. 

Case 2: There exists an angle in triangle ABC greater than or equal to 120°. 
Only one angle of the triangle is greater than or equal to 120°, say ZACB > 120°, see Figure [2jb) for 
example. Then the Fermat point is coincident with C. 

4 The MMA problem 

We solve the MMA problem depending on the degree of vertex v. 
4.1 Vertex of Degree 2 

Suppose that the degree of v is two. Let A and B be the positions of its adjacent vertices. The task 
is to find a point P* with |PP*| < r maximizing angle 9 = mm(ZAP'B,ZBP'A) (the angles are in the 



3 



counterclockwise order). Suppose that segment AB intersects circle Z. Obviously, 8 = 180° is achieved by 
any point P* in the intersection of AB and Z. The solution is unique if ABHZ is a single point (for example, 
AB may be tangent to Z). 




(a) (b) 

Figure 3: The case where AB n Z ^ 0. (a) P* = Q where Q is defined by Equation ([!]). (b) Point Q lies 
not in AB n Z. 

Suppose that ABdZ is a segment of positive length, see Figure [3]( a). There are infinitely many positions 
of P* (all in AB n Z achieving 6 = ir). If P is not on segment AB then any location of P* will change 
distances to A and B. We would like to preserve the ratio |AP|/|I?P| and compute point Q in segment AB 
such that \AQ\/\BQ\ = \AP\/\BP\. This implies 

\AQ\ ^ \AP\ 

\AB\ \AP\ + \BP\' [ ' 

Then the coordinates of point Q can be computed using 



Qx — A x 



Qy ~ Ay 



\AP\ 
\AP\ + \BP\ 

\AP\ 
\AP\ + \BP\ 



{B x -A x ), 



(By~ A y ). 



(2) 
(3) 



If point Q lies in circle Z then P* = Q as shown in Figure ^a) . Otherwise we select the endpoint of 
segment AB n Z that is closer to Q, see Figure [3^b) for an example. 
We assume now that AB and Z do not intersect. 

Proposition 1 If segment AB and circle Z do not intersect then P* is on the boundary of Z. Furthermore 
the circle passing through A,B, and P* is tangent to Z. 

Proof: Suppose P* — (x,y) where x and y are unknown. Consider circular arc AP*B as shown in Figure 
|4|a). The angles ZAP' B are equal for all points P' from arc AP* B (since the inscribed angle of a chord 
equals half of the central angle ,see Figure [§»). Therefore we can assume that P* is on the boundary of 
Z. Furthermore arc AP*B cannot intersect the boundary of Z twice since any point P' from the arc of Z 
cut off by arc AP* B makes an angle larger than P* makes, i.e. ZAP' B > ZAP*B as shown in Figure [4] 
Thus, arc AP* B must be tangent to Z as shown in Figure [5ja). ■ 

The MMA problem can be formulated now as the following problem dealing with only one angle. 

Maximum Angle Problem. Let A. B, and P be points in the plane and a number r > such 
that \AP\, \BP\ > r. Compute a point P* with |PP*| = r maximizing angle ZAP*B. 

Let O be the center of the circle passing through A, B, and P* . It suffices to find point O since point P* 
is the intersection point of the circle Z and the line passing through points O and P. The first equation is 
(P* — P x ) 2 + (P* - P y ) 2 = r 2 and the second equation is linear 
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where the subscripts denote the coordinates. 

Using coordinate transformations we can assume without loss of generality that ^4(0,-1) and -8(0, 1), 
see Figure |5jb). Then \OA\ = \OB\ implies that O(x,0) where x is unknown. Let r' = \OA\ be the radius 
of the circle centered at O. Then \OP\ = r + r' and \OA\ = r' . They can be written as 

U x -P x f + pi = {r + r'f (4) 
\ x 2 + l= r ' 2 (5) 

Subtracting Equation Q from Equation ([5]) , we obtain 

2P x x = Pi + Py -2/r-r 2 - 1. (6) 

By plugging x from this equation into Equation of ([5| we obtain a quadratic equation in the variable 
r' . There are two roots of the quadratic equation and they correspond to two circles shown in Figure [6] 
Therefore we take the smallest root of the quadratic equation. 

Therefore we proved the folowing claim. 

Proposition 2 The maximum angle problem can be solved using a quadratic equation. 



5 Vertex of degree 3 

In this Section we consider the case where the degree of v is three. Let A, B and C be the positions of its 
adjacent vertices in G. We consider the first case of the Fermat point of triangle ABC where all angles of 
ABC are less than 120°. 

Proposition 3 // all angles of ABC are less than 120° and the Fermat point F of ABC lies within circle 
Z, then P* = F. 
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Figure 6: An example of the maximum angle problem. The circles centered at O and O' are the two solutions 
of the system of equations Q and ^ . 



Proof: Without loss of generality we assume that A,B,C are in counterclockwise order around P* , see 
Figure [7] (a). Suppose to the contrary that P* is a point inside Z different from F as shown in Figure [7] 

(b) . The Fermat point lies in one of triangles AP* B, BP*C, or AP*C. Suppose that it belongs to triangle 
AP*C. Then ZAP*C < ZAFC = 120° and P* is not the optimal point. Therefore P* = F, see Figure [7] 

(c) .- 




(a) (b) (c) 

Figure 7: (a) Setting of Proposition [3] (b) Contradiction when P* ^ F. (c) Points P* and F coincide. 

Note that in the case of Proposition [3] all angles around P* are equal 120°. In the rest of this section 
we consider cases where not all angles ZAP* B, ZBP*C, and ZAP*C are equal. If the smallest angle of 
ZAP*B, ZBP*C, and ZAP*C is unique, say it is ZAP*B, then it can be computed by solving the maximum 
angle problem for A, B using the method from the previous section. It remains to consider the case of two 
smallest angles. 



5.1 Two smallest angles 

In this section we assume that there are exactly two smallest angles from {ZAP*B, ZAP*C, ZBP*C}. First, 
we show that P* must be on the boundary of circle Z. 

Lemma 4 If the degree of v is 3 and there are two smallest angles around P* then P* lies on the boundary 
ofZ. 

Proof: We assume that 

(i) A,B,C are in counterclockwise order around P* , and 

(ii) ZAP*B = ZBP*C < ZCP*A. 

Suppose to the contrary that P* lies inside Z as shown in Figure [8] Draw two circles, one passing through 
points A,B,P* and the second passing through points B,C,P*. Since these circles intersect by two points 



G 



B and P* , then their interiors intersect by a lune. Since P* is inside circle Z, then the intersection of the 
lune and the interior of Z is a non-empty set /, see the shaded area in Figure |sja) . For any point P' in /, 
ZAP'B > ZAP*B and ZBP'C > ZBP*C (this can be seen, for example, by the fact that the inscribed 
angle of a chord equals half of the central angle, see Figure |8jb)). We choose P' from I close enough to P* 
such that ZAP'B, ZBP'C < ZCP'A. Then the smallest angle around P' is greater than ZAP*B. Therefore 
P* is not the solution of the MMA problem. Contradiction. ■ 




(a) (b) 

Figure 8: (a) For all points P' in the shaded area ZAP'B > ZAP* B and ZBP'C > ZBP*C. (b) The 
inscribed angle ZAP*B of chord AB equals half of the central angle ZAOB. 

The main result in this section is the following 

Theorem 5 If the degree of vertex v is three and there are two smallest angles around P* then P* can be 
computed by solving a polynomial equations of degree at most four. 




Figure 9: Points Pi and Pi are the points on circle Z maximizing angles ZAP\B and ZBP2C, respectively. 
Point P* is different from P x and P 2 and angles ZAP* B = ZBP*C. 

Proof: First we show that it is possible that there are two smallest angles. An example where two smallest 
angles around P* are equal is shown in Figure [9] where points P, P* and B are collinear and points A and 
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C are symmetric about the line PB. Therefore angles ZAP 2 B < ZAP* B < ZAP 4 B and P* is the solution 
of the MMA problem. 

We now prove the theorem. Suppose that ZAP*B = ZBP*C < ZCP*A. We need to find coordinates 
of P*(x,y). By the law of cosines we write the equation cos(ZAP* B) — cos(ZBP*C) as 

\BP*\ 2 + \AP*\ 2 - \AB\ 2 _ \BP*\ 2 + \CP*\ 2 - \BC\ 2 

2 • | BP* | • \ AP*\ ~ 2 • | BP* | • | CP* | *■ ' 

\CP*\ ■ (\BP*\ 2 + \AP*\ 2 - \AB\ 2 ) = \AP*\ ■ (\BP*\ 2 + \CP*\ 2 - \BC\ 2 ) (8) 

Without loss of generality we can assume that P = (0,0) and A y — C y . By Lemma H we assume 
|PP*| = r. Then x 2 + y 2 = r 2 and 

|CP*| 2 = r 2 + C 2 + C 2 - 2C x x - 2C v y = a x x + a 2 y + a 3 (9) 
X - (|BP* | 2 + |AP*| 2 - | AB\ 2 ) =r 2 + A X B X + A y B y - (A x + B x )x - {A y + B y )y = a 4 x + a 5 y + a 6 (10) 



\AP* \ 2 =r 2 +A 2 x + A y - 2A x x - 2A y y = b x x + b 2 y + b 3 (11) 
(|BP*| 2 + |P*C| 2 - |BC| 2 ) = r 2 + B X C X + B y C y - (B x + C x )x - (B y + C v )y = b 4 x + b 5 y + b 6 , (12) 



where <n = -2C x ,a 2 = -2C y ,a 3 = r 2 + C 2 + C y 2 ,a 4 = -(A x + B x ),a 5 = -{A y + B y ), 

a 6 = B X A X + B y Ay + r 2 , h = -2Ax, b 2 = -2A y ,b 3 = r 2 + A 2 + Ay 2 ,b 4 = -{B x + C x ),b 5 = -{By + Cy), 

b 6 = B X C X + B y C y + r 2 . 

Take the square of both sides of Equation ^ 

(aix + a 2 y + a 3 )(a 4 x + a 5 y + a 6 ) 2 = (biX + b 2 y + b 3 )(b 4 x + b 5 y + b 6 ) 2 (13) 

It can be written as 

cix 3 + c 2 y 3 + c 3 x 2 + c 4 y 2 + c 5 x 2 y + c 6 xy 2 + c 7 x + c$y + c^xy + c w = (14) 

where c\ = a\a\ - b\b\ c 2 = a 2 a\ - b 2 b\, c 3 = (2a 4 a 4 a e + a 3 a| - (2bib 4 b e + b 3 b\)), 
c 4 = (2a 2 a 5 a 6 + a 3 a 2 - (2b 2 b 5 b 6 + b 3 b\)),c b = (2aia4a 5 + a 2 a 2 - (2bib 4 b 5 + b 2 b 2 4 )), 
c 6 = (2a 2 a 4 a 5 + a 4 a 2 — (2b 2 b 4 b 5 + b\b\)), c 7 = (2a 3 a 4 a 6 + a x a\ — (2b 3 b 4 b 6 + b\b%)), 

c 8 = (2a 3 a 5 a 6 + a 2 a\ - (2b 3 b 5 b 6 + 6 2 &1))j c 9 = 2(a 1 a 5 a 6 + a 2 a 4 a 6 + a 3 a 4 a 5 - (&1&5&6 + b 2 b 4 b 6 + b 3 b 4 b 5 )), 
Cio = a 3 ag. 

Using y 2 = r 2 — x 2 we reduce powers of y in Equation ( 14 ) 



Cix 3 + y{c 2 r 2 + c 8 ) + c 3 x 2 + c 4 — c 4 x 2 + x 2 y(c 5 — c 2 ) + c 6 x — c 6 x 3 + c 7 x + c s y + c 9 xy + c 10 = (15) 

d!X 3 + d 2 x 2 + d 3 x + d 4 + (d 5 x 2 + d 6 x + d 7 )y = (16) 

where d\ = ci — c 6 , d 2 — c 3 - c 4 , d 3 = c 6 r 2 + c 7 , d 4 = c 4 r 2 + c W: d 5 = c 5 - c 2 , ^6 = eg, d 7 = c 2 r 2 + c 8 . 



(dix 3 + d 2 x 2 + d 3 x + d 4 ) - ((d 5 x 2 + d 6 x + d 7 ) 2 ■ (r 2 - x 2 )) = (17) 
Simplifying Equation ( fl7| we get a sextic equation 

eix 6 + e 2 x 5 + e 3 x 4 + e 4 x 3 + e^x 2 + e§x + e 7 = (18) 

where e x =d\+ d 2 5 ,e 2 = 2{d x d 2 + d 5 d e ), e 3 = {d 2 2 + 2d x d 3 - (d 2 5 r 2 - (d 2 6 + 2d 5 d 7 ))), 

e 4 = (2(d 1 d 4 + d 2 d 3 ) - {2(d 5 d 6 r 2 - d 6 d 7 ))),e 5 = ((d 2 + 2d 2 d 4 ) - {(d 2 + 2d b d 7 )r 2 - d 2 7 )), 

e 6 = 2d 3 d 4 — (2d e d 7 r 2 ),e 7 = d\ — d 2 r 2 . 

By Lemma [6| the sextic equation has a quadratic factor. Dividing by it, we obtain a polynomial of degree 
four (Equation (|19|,(20), or (21)) and the theorem follows. ■ 
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Lemma 6 The polynomial equation (18) has factor x 2 + A 2 — • 



Proof: We consider 3 cases. 



Case I: A y = r. Then es = e§ = e-j = and the polynomial ( 18 ) has factor x 2 . The polynomial is reduced 



to the polynomial of degree 4 



e\x A + e 2 x 3 + e 3 x 2 + e±x = 0. 



(19) 



Case II: A 2 — r 2 > 0. We scale the coordinates such that A 2 — r 2 — 1. Then e 6 



e 4 - e 2 , e 7 = e 3 - e x 



and the polynomial (18 1 has factor x + 1. The polynomial is reduced to the polynomial of degree 4 
e\X A + e 2 x 3 + (e 3 - ei)x 2 + (e 4 - e 2 )x + e 5 - (e 3 - ei) = 0. 



(20) 



Case III: A 2 ~r 2 < 0. We scale the coordinates such that A 2 — r 2 = —1. Then cq — e 4 + e 2 , ej = e 3 + e\ 
and the polynomial (18 1 has factor x 2 — 1. The polynomial is reduced to the polynomial of degree 4 



e\X A + e 2 x 3 + (e 3 + ei)x 2 + (e 4 + e 2 )x + e 5 + (e 3 + ei) = 0. 



(21) 



The lemma follows. 



6 Algorithm 

In this Section we discuss how to modify the spring embedder in order to take into account the angles between 
embedded edges. Let G = (V, E) be an input graph. At the initial step the spring embedder randomly places 
the vertices of G in the plane. Then, it iterates a simultaneous motion of the vertex positions based on one 
or several forces (springs). We describe a new force using angle optimization. For every vertex v with a 
position P in the plane, the algorithm computes a new position P* and uses vector PP* as a force applied 
to vertex v. 

Angle Optimization Algorithm 
Input: Graph G = {V,E) embedded in the plane and radius r. 

Output: New embedding of graph G where each vertex is translated within distance r. 
For each vertex v £ V do the following steps. 

1. Let P be the position of v in the plane. We compute P* as follows. Every time P* is assigned, we 
proceed to the next vertex v. 

2. If the degree of v is at most one then set P* = P. 

3. Suppose that the degree of v is equal to two. Let A and B be the positions of vertices adjacent to v. 

(a) Suppose AB n Z ^ 0. Compute point Q using Equations ^ and If \PQ\ < r then P* = Z: 
otherwise assign P* to the endpoint of segment AB n Q that is closer to Q. 

(b) If AB n Z = then set compute P* using Proposition [5] 

4. Suppose that the degree of v is equal to three. Let A,B, and C be the positions of vertices adjacent 
to v. Compute Fermat point F of triangle ABC. If \PF\ < r then set P* = F. Otherwise, for each 



segment ab £ {AB, AC, BC}, compute point P a b maximizing angle ZaP^b (Section 4.1). If angle 
ZaPabb is the smallest angle from {ZAPabB, ZAPacC, ZBPbc'G} then set P* = P a \,. Otherwise 
compute P* as the best solution using two smallest angles for every two segments from {AB, AC, BC} 



(Section 5.1 ) 



For the remaining vertices of degree at least four, apply the following grid approach. Pick a grid 
stepsize 6, for example S = r/3. Consider a grid with the origin at P. For every grid point Q with 
\PQ\ < r, compute the smallest angle q.q if P is moved to Q. Find Q maximizing qq. Assign P* = Q. 
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(a) 



(b) 



Figure 10: Graph Tiq from [17] . (a) Drawing by Spring Embedder pQ. (b) Drawing by our program with 
angle optimization. 



We implemented this algorithm^] and run it on several graphs. First, we tested the algorithm on graph 
T 10 from |17j since it contains vertices of degree two. The program draws T 10 with angles optimized, see 
Figure [To] (b) . It can be compared with the drawing of the original embedder [T] , see Figure [To] (a) . 

We also tested the algorithm on the phylogenetic networks for Algae example [IBJ. The drawing of the 
Algae network by the spring embedder pQ is shown in Figure [I] It can be compared with our drawings in 



Figure 11 In all drawings (in Figures!!] and Figure 11) the edge length constraints are satisfied but the 



angle resolution in drawings in Figure [Tl| is significantly larger 



shown in Figure 11 (a) uses the weighted version of the graph and the drawing in Figure 11 (b) uses the 
graph with intermediate points on the edges. 





(a) 



(b) 



Figure 11: Phylogenetic network for the Algae example [TB]- (a) Weighted graph, (b) Graph with interme- 
diate points on the edges. 



Finally, we run our program on the well-known graphs: the Petersen graph [13 , the Heawood graph [12j . 



and the Herschel graph [5], see Figure 12 The Petersen graph is drawn with two crossings (the Petersen 
graph is in fact the smallest 2-crossing cubic graph), see Figure [l2^a). The Heawood graph has crossing 
number three (it is actually the smallest 3-crossing cubic graph) and our program found a drawing with 



exactly three crossings, see Figure 12 'b) 



1 Demo is available at http://www.utdallas.edu/~sxb027100/soft/Angle0pt/ 
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(a) (b) (c) 

Figure 12: (a) The Petersen graph, (b) The Heawood graph, (c) The Herschel graph. 

7 Conclusion 

We proposed a novel approach to the problem of optimizing the angular resolution of a drawing. It has been 
applied to the spring embedder and the results with good angular resolution were obtained. It is known that 
the running time of the spring embedder can increase with the size of the graph. Therefore it is important to 
perform better the intermediate steps. The optimal solution of the MMA problem provides an opportunity 
to decrease the number of iterations of the spring embedder. 

The main result of this paper states that a vertex of degree at most three can be displaced optimally by 
solving a polynomial equation of degree at most four (which is interesting since the straightforward approach 
leads to a polynomial of degree six) . A special case of MMA problem for degree three verices is the classical 
Fermat problem and the Fermat point is the solution for the special case. An interesting question for future 
research is to find the algebraic complexity of MMA problem for higher vertex degrees. 
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